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ABSTRACT 

We present an analysis of the neutron star High Mass X-ray Binary (HMXB) 4U 1909+07 mainly based on Suzaku data. We extend 
the pulse period evolution, which behaves in a random-walk like manner, indicative of direct wind accretion. Studying the spectral 
properties of 4U 1909+07 between 0.5 to 90 keV we find that a power-law with an exponential cutoff can describe the data well, when 
additionally allowing for a blackbody or a partially covering absorber at low energies. 

We find no evidence for a cyclotron resonant scattering feature (CRSF), a feature seen in many other neutron star HMXBs sources. 
By performing pulse phase resolved spectroscopy we investigate the origin of the strong energy dependence of the pulse profile, 
which evolves from a broad two-peak profile at low energies to a profile with a single, narrow peak at energies above 20keV. Our 
data show that it is very likely that a higher folding energy in the high energy peak is responsible for this behavior. This in turn leads 
to the assumption that we observe the two magnetic poles and their respective accretion columns at different phases, and that these 
accretions column have slightly different physical conditions. 
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1. Introduction 

High mass X-ray binaries (HMXBs), in which the compact ob- 
ject is a rotating neutron star, provide a unique way to investi- 
gate the flow of matter into a deep gravitational potential. By 
analyzing the evolution of the pulse period and its correlation to 
luminosity, the presence or absence of a stable accretion disk can 
be probed (Ghosh & Lamb 1979| . By performing pulse phase 
resolved spectroscopy, different areas of the X-ray producing re- 
gion can be analyzed and thereby the interaction between ra- 
diation and matter in very strong magnetic fields can be stud- 
ied. 4U 1909+07 is an ideal source for these investigations, as 
it is a slowly rotating, persistent HMXB. Discovered by Uhuru 
( |Giacconi et al. 1974), its 4.4 d o rbital period was determined in 
RXTE/ASM data by |Wen et al.| ([2000). Pulsations with a pulse 
period of «605 s were later found in RXTE/PCA data by|Levine 



|et al. (2004 1. These authors also found that the intrinsic absorp- 
tion Nh is strongly variable over the orbit, peaking around orbital 
phase 0. These measurements can be explained by a smooth stel- 
lar wind and an inclination of the system between 38° and 68° 
( |Levine et alj |2004| >. |Fiirst et al.| ( |2011a| hereafter Fll) investi- 
gated the long-term evolution of the pulse period between 2003 
and 2009 using INTEGRAL/lSGRl lightcurves. We found that 
the evolution is consistent with a random walk, a typical behav- 
ior for a wind accreting source without a stable accretion disk, 
similar to Vela X-l (see, e.g., de Kool & Anzer T993J. 



The spectrum of 4U 1909+07 can be described using typical 
phenomenological models often applied to HMXB neutron star 
sources, such as a power-law with an exponential cutoff (Levine 
|et al.| [2004] Fll). By analyzing phase-resolved RXTE data we 



found that a blackbody component with a temperature around 
1.7keV is also needed to describe the soft excess in the spec- 
trum. Many neutron star sources show a cyclotron resonant scat- 
tering feature (CRSF, see, e.g., |Schonherr et al. 2007|> at en- 



ergies between 10 and 50keV. Neither in RXTEfHEXTE nor 
in iNTEGRAL/ISGRI data is evidence for such a feature seen; 
however, the signal-to-noise ratio (S/N) in these data is rather 
low. To investigate the behavior of the soft excess in more de- 
tail and to perform a rigorous search for a CRSF we obtained a 
20 ks Suzaku observation, the analysis of which is presented in 
this article. In Sect. [2] we give a short overview of the data used 
and the reduction pipeline. In Sect. [3] we discuss the evolution of 
the pulse period, while in Sect.|4]we analyze the X-ray spectrum. 
We perform phase averaged (Sect. 
spectroscopy (Sect. 4.2 1. In Sect, 
our results. 



4.1| i as well as phase resolved 
51 we summarize and discuss 



2. Observations & data reduction 

This paper is mainly focused on a 20 ks Suzaku observation per- 
formed on November 2, 2010. Data were taken in the orbital 
phase range 0.472-0.627, i.e., outside the phases of increased 
absorption. We use data from the X-ray Imaging Spectrometers 
|2007| l 0, 1, and 3 in the energy range 0.5 



(XISs, Koyama et al. 



lOkeV and from the High X-ray Detector's (HXD, Takahashi 



|et~aTl [2007) 1 PIN instrument from lOkeV to 90keV. We did not 
use data from the HXD/GSO due to source confusion with the 
nearby, bright black-hole binary GRS 1915+105. As the hard 
X-ray spectrum of 4U 1909+07 falls off exponentially, a sig- 
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nificant detection in GSO would be rather unlikely. The data 
were reduced using the standard pipeline from HEASOFT v6.1 1 
and the HXD calibration from September 201 1, including dead- 
time correction for PIN. Data analysis was performed using 
the Interactive Spectral Interpretation System (ISIS, |Houck &| 
Denicola) |2000| l, version 1.6.2. We extracted lightcurves from 
XIS with 16 s time resolution in the energy bands 0.5-5 keV, 
5-10 keV, and 0.5-10 keV, as well as from PIN with 40 s res- 
olution in the energy bands 10-20keV, 20-40 keV, and 40- 
90keV. Figure [T] shows the lightcurves from all instruments, re- 
binned to 160 s for better visibility. Some flaring is evident in the 
lightcurve, however, no significant changes in hardness could be 
found when comparing the different energy bands. For compar- 
ison we also extracted data from the ISGRI detector ([Lebrun 
eTaT|[2003l ) aboard INTEGRAL ( |Winkler et al] [20031 ), between 
MJD 55107.13 to 55148.82 and MJD 55503.74 to 55518.02 
with exposure times of 235 ks and 396 ks, respectively. Using 
the standard OSA 9.0 we extracted lightcurves between 20 to 
40keV with 20 s time resolution. All lightcurves, Suzaku and 
INTEGRAL, were transferred to the barycenter of the solar sys- 
tem and corrected for the orbital motion of 4U 1909+07, using 
the ephemeris by Levine et al. (2004|. 



3. Timing 

The pulse period evolution of 4U 1909+07 was published by 
Fll up to spring 2009. Figure [2] extends this time series up to 
fall 2010, using all publicly available INTEGRAL data and the 
Suzaku observation. INTEGRAL also observed 4U 1909+07 
in spring 2010, but the rough sampling does not allow for 
a reliable pulse period measurement. To obtain the pulse pe- 
riods we performed epoch folding ( |Leahy et al.| |1983[ ) us- 
ing XIS lightcurves between 0.5-10keV with 16 s time resolu- 
tion, PIN lightcurves between 10-90 keV with 40 s time reso- 
lution, and INTEGRAL /ISGRI lightcurves between 20^10 keV 
with 20 s time resolution. Uncertainties on the pulse periods 
of INTEGRAL were obtained by fitting a Gaussian distribu- 
tion to the x 1 distribution of the individual epoch folding re- 
sults, as described by Fll. The epoch folding distribution of 
Suzaku/XIS as well as Suzaku/PIN does not show a well enough 
defined peak for fitting with a Gaussian function. We therefore 
estimated the uncertainties on the period measurement using a 
Monte Carlo approach as described by Davies ( 1990| l, where 



we simulated 1000 lightcurves with the same binning as the 
original measurement, scattered around the measured pulse pro- 
file. The width of the distribution of resulting pulse periods in 
these lightcurves is a good indicator of the uncertainty of the 
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Fig. 2. Pulse period evolution between 2003-2011. Black data 
were already published by Fl 1 . More recent INTEGRAL /ISGRI 
data are shown in red, Suzaku/PIN in blue. Suzaku/XIS is not 
shown to ensure readability. 



pulse period. This approach leads to Pxis = 604.10 ± 0.20 s and 
fpiN = 603.87 + 0.24 s (all uncertainties are given at the 90% 
level). Due to the longer exposure time, the almost simultane- 
ous INTEGRAL /ISGRI data provide a much better constrained 
pulse period measurement, with Pisgri = 603.86±0.05 s. As this 
period is consistent with the Suzaku measurements, we will use 
it throughout this paper and for the analysis of the Suzaku data. 

The pulse profiles of Suzaku/PIN and INTEGRAL /ISGRI 
in the 20^1-0 ke V range agree, showing a short single peak fol- 
lowed by a long dim phase (Fig. [3}l). This profile is in agree- 
ment with RXTE/HEXTE (see, e.g., |Levine et al."t[2004] Fll). 
As described by Fl 1 a strong change in the pulse profile occurs 
between 8 and 20 ke V, where it changes from a two-peaked pro- 
file to a single peaked profile. This energy dependent behavior is 
also seen in other neutron star X-ray binaries, like 4U 01 15+634 
( |Mulleretal.|[20T0| and 1A 1118-61 ( |Suchy et al.|[20TT) . The 
energy range of the transition is only marginally covered by the 
instruments aboard Suzaku, as it is just in the overlapping re- 
gion between XIS and PIN. The XIS pulse profile in different 
energy bands is shown in Fig. [3^ and Fig. [3|? and is very similar 
to RXTE/PCA profiles as presented by |Levine et al.| ( |2004| l and 
Fll. 



4. Spectroscopy 

4. 1 . Phase-averaged spectroscopy 

As discussed by Fll, the phase averaged spectrum of 
4U 1909+07 can be very well described with typical phe- 
nomenological models often used to describe the spectra of neu- 
tron stars in HMXBs plus a Gaussian-shaped fluorescent Fe Kcc 
line. T he XSPEC m odels cutoffpl, highec ut ( |White et aT) 
|1983| ), and comptt ( |Hua & Titarchuk[ |1995| > all gave similar 
statistically acceptable fits. Therefore, we applied these mod- 
els to the Suzaku data and the results are shown in Tables [T] [2] 
and|3] respectively. The data are shown together with the best-fit 
cutoffpl model in Fig. |4] 

The flux of the Suzaku data (XIS normalizatiorQ is compa- 
rable to the RXTE data (PCA normalization) presented by Fl 1, 
with^ s ° u ™ eV = 0.324™ keVs- 1 cor 2 and ^ T » )0keV = 
0.330 ± 0.007 keV s _1 cirT 2 , respectively. The Fe Ka line is nar- 
row and its width could not be resolved in the XIS data. We 



PIN cross-calibration factor was within 2% 
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Fig. 3. Energy resolved pulse profiles as measured with 
Suzaku/XIS (a) and b)) and PIN (cj to ej). Marked in color and 
labeled A-E are the phase bins used for phase-resolved spec- 
troscopy (see text for details). In panel c), INTEGRAL /ISGRI 
data is superimposed in red, scaled according the right-hand y- 
axis. 
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Fig. 4. a) XIS0 (red), XIS1 (green), XIS3 (orange), and PIN 
(blue) spectrum together with the best fit cutoffpl-model in- 
cluding partial covering. In magenta the strongly absorbed part 
is shown, in dark blue the less absorbed one. The Fe Ka line is 
shown in brown, b) residuals to the cutoffpl model without 
partial coverer. c) residuals to the best-fit model. 



Table 1. Phase averaged spectral parameters for the cutoffpl 
continuum model. In the second column the best-fit parameters 
with an additional blackbody component (bbody) are shown, in 
the third the parameters when including a partial coverer (PC). 
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therefore fixed it to 10~ 6 keV in all fits. It was shown by Fll 
that an additional blackbody component is necessary to describe 
a soft excess. Adding this component leads also to a clear im- 
provement in the fit to the Suzaku data, as can be seen in the 
second column of Tables [T] |2j and[3] However, some continuum 
parameters become only very weakly constrained, like Ef \d in 
the cutoffpl and highecut model or the seed photon tem- 
perature r in the comptt model. Additionally the 20-40 keV 
flux drops to unrealistically low values^] and the best fit black 
body temperature is very different between the cutoffpl and 
highecut model, rendering a physical interpretation impossi- 
ble. To describe the soft energy part of the spectrum without 
a black body component, we included a partial coverer in the 
model with two different absorption columns /V H ,i an d Nu,2^ with 
their relative influence determined by a covering fraction CF: 

(CF x 7%i + (1 - CF) x /V H ,2) x (cont + Fe Ka) (1) 

This model leads to equally large improvements in terms of^- 2 , 
as seen in the third column of Tables [T] [2] and [3] Although 
the secondary absorption column Nh,2 is also only weakly con- 
strained, the continuum parameters remain very well constrained 
and at physically reasonable values. 



4.2. Phase-resolved spectroscopy 

To study the variance of the spectrum with the viewing angle 
onto the neutron star, i.e., the pulse phase, we divided the data 
in 5 phase bins, as indicated in Fig. [3] The bins were chosen 
to cover the prominent features, e.g., the deep minimum be- 
tween the pulses and the hard primary pulse (bins C and D, 
respectively) as well as to retain a S/N sufficient to constrain 
the individual spectral parameters well. As shown by Fll, the 
cutoffpl model delivers the best description of the phase- 
resolved spectra, and in order to be able to compare our results 
directly with the previous work, we applied that model to the 
phase-resolved data. Neither a blackbody component nor a par- 
tial coverer could be significantly detected in these spectra due to 
the lower S/N, therefore we only use a single absorber, the con- 
tinuum, and a Gaussian-shaped fluorescent iron line at 6.4 keV. 



compensated by a very high PIN cross-calibration factor of asl.3 
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Table 2. Same as Tab.[l] but for the highecut model. 
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Table 3. Same as Tab.[l] but for the compTT model. 



Detailed analysis of the phase averaged spectrum shows that 
there are strong systematic correlations between YVh, F, and Ef \&. 
These correlations lead to large uncertainties in the phase re- 
solved spectra in all parameters. To reduce these uncertainties 
we fixed each of the three parameters in turn, and compared the 
results with the energy resolved pulse profiles. All fits lead to 
very similar ^ ed values, so that a distinction based on the quality 
of the fit is not possible. However, when comparing the fluxes 
in the 2CM-0 keV energy band modeled with a fixed folding en- 
ergy to the pulse profile in that energy range, a clear discrep- 
ancy is evident, i.e., the bright main peak is not correctly re- 
produced. Additionally, Fl 1 found indications that the phase re- 
solved RXTE spectra are best described using a variable Ef n. 
We therefore exclude the model with a fixed folding energy from 
the following discussion. 

As can be seen in Fig. [5J, the folding energy is in all mod- 
els very high during phase bins C and D, the minimum and the 
primary peak. The folding energy seems therefore responsible 
for the spectral hardening during the primary peak, and not, e.g., 
strong changes in the photon index T. The results also show that 
the minimum has a similar spectrum as the primary peak. In the 
models with variable T, stays high during phase bin E, but 
at the same the time T is rising, leading to an overall softer spec- 
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Fig. 5. Results of the phase resolved spectroscopy. In black the 
fit without any fixed parameters is shown, for the blue data points 
the Nh was frozen, for the red the photon index T. In green 
the results of the phase averaged fit are shown, a) XIS0 pulse- 
profile in the 0.5-10 keV range, b) flux between 3-10 keV in 
keVs~'cirr 2 , c) flux between 20^-0 keV in keVs _1 cirr 2 , d) 
absorption column density, e) photon index, f) folding energy, 
g) flux of the iron line in ph s _I crrT 2 . 



trum. With F frozen, Ef^ drops to values similar to those in 
phase bins A and B. Differences between these two descriptions 
would only be visible above ~70 ke V, where the source flux has 
dropped below the detection limit of PIN. 

Variations of YVh are visible in Fig. [5}i, but seem uncorre- 
lated with the observed hard primary peak. YVh drops to low val- 
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ues during the minimum in phase bin C, but the overall variation 
is only marginally significant. It is interesting to note that the 
flux of Fe Ka is also variable with phase (Fig.[5j*), however it is 
not correlated with the overall soft X-ray flux nor the A^h varia- 
tions. There rather seems to be a slight shift with the Fe Ka flux 
reacting on the 3-10keV flux with approximately 120s delay. 
Using the model in which all parameters were allowed to vary 
and shifting the Fe Ka flux by one phase bin, the correlation 
to the soft continuum increases from 0.04 to 0.70 according to 
Pearson's correlation coefficient. The energy of the Fe Ka line 
(not shown) is not varying significantly with pulse phase. 

5. Summary & Outlook 

We have presented a timing and spectral analysis of a 20 ks 
Suzaku observation of 4U 1909+07. We extended the pulse pe- 
riod evolution as presented by Fll by two data points, the 
most recent of which, P — 604.86 s, is consistent between 
INTEGRAL and Suzaku. To describe the spectrum we applied 
different phenomenological models and find statistically accept- 
able fits when including a fluorescent Fe Ka line at 6.4 keV and 
either a blackbody component or a partially covering absorber. 
The continuum parameters in the models with a blackbody com- 
ponent are very similar to the ones found by Fll in RXTE- 
data. However, they are highly unconstrained in the Suzaku- 
data, while using a partial coverer still results in a good descrip- 
tion of the continuum with small uncertainties. 

A partial coverer can be used as an approximation to a com- 
plex structured stellar wind ("clumpy" wind), as seen in many 
HMXBs as well as isolated O and B stars (see, e.g., Oskinova] 



et al. 2007 Fiirst et al. 201 lb| among many others). Within 
these winds, dense structures with column densities of the order 
of 10 24 cirT 2 can be reached, as seen in many other sources (e.g., 
GX 301-2 o r IGR J 163 18-4848; |Fiirst et al.| |2011b[ |Barragan| 
et al.| |2009[ l. The absorption column of the partial coverer in 
4U 1909+07 is measured to be also of that order, with values 
between (0.5— l)xl0 24 cirT 2 , as seen in Tables [I] [5] and|] 

We also applied the partial covering model to the older 
RXTE data, which have about twice the column density as the 
Suzaku data, but did not find an acceptable fit. As the RXTE 
data do not cover the very soft energies, a small contribution of 
the partial coverer can easily be missed in the data. If the partial 
coverer is due to the stellar wind, it is also highly variable, as 
structures will move with th e wind with typical outflow termi- 
nal velocities > 1000km s _1 (Prinja et al. 1990 1. Therefore, the 
RXTE observation might have just been performed during a less 
obstructed view onto the X-ray source. Nonetheless, these data 
clearly require a blackbody component to describe the data. 

Assuming that the model including the partial coverer repre- 
sents a good description of the physical conditions in and around 
the X-ray source, it means that the blackbody disappeared, i.e., 
that it is highly variable. The blackbody component originates 
from the thermal mound on the neutron star or the accretion col- 
umn and its properties can vary with changes in the accretion 
rate and the physical conditions in the accretion column close 
to the neutron star. The X-ray fluxes of the Suzaku and RXTE 
data are, however, similar, i.e., it is likely that the accretion rate 
and thereby the temperature in the column are also very simi- 
lar. Changes in the geometry of the accretion column could ac- 
count for small differences, like the ones seen when comparing 
the RXTE model to the best fit Suzaku model which includes a 
blackbody component. In Suzaku the temperature is measured to 
be around 2 keV, compared to « 1.5 keV in the RXTE data. That 
the blackbody component becomes completely invisible due to 



changes in the geometry of the accretion column seems unlikely. 
It is also possible that the blackbody is a spurious effect in the 
RXTE data. For example, the galactic ridge emission, which can 
strongly influence RXTE/PCA data due to the large field of view 
of the collimator, can also often be described by a soft thermal 
component (see, e.g., Yamauchi et al.| 2009 [ l. The galactic ridge 
component is suppressed in Suzaku/XIS due to its imaging capa- 
bilities. To study the variability of the blackbody and the partial 
coverer, a more extensive monitoring in the soft X-rays is neces- 
sary. 

To investigate the behavior of the spectrum under different 
viewing angles, we performed phase resolved spectroscopy in 5 
phase bins. As already seen in the energy resolved pulse profiles 
(Fig. [3}, the primary peak (labeled D) is distinctly harder than the 
broad secondary maximum (bins A, B, and E). We found that 
this is likely due to a reduced folding energy in the secondary 
maximum. This change can be explained by assuming that ac- 
cretion happens on both poles of a magnetic dipole field and 
that each pulse maximum corresponds to one of the poles. Small 
differences in magnetic field strengths will lead to different ge- 
ometries and thereby temperatures in the two accretion columns. 
The hotter the electron plasma the more effectively it can scatter 
photons to high energies via the inverse Compton effect, leading 
to an overall harder spectrum. This picture might be over sim- 
plified, as relativistic light bending effects and the geometry of 
the accretion column need to be taken into account, but such an 
analysis is beyond the scope of this article. There are, however, 
new models under development to describe the X-ray producing 
region more realistically (e.g., |Becker et aL 2012 1, which we 
will apply in a future publication. 

We searched for CRSFs in the phase averaged as well as in 
the the phase resolved spectrum, but did not find any evidence 
for such a feature in the energy range between 10 and 70keV. 
Assuming typical values of £crsf = 25keV and cr CRSF = 
10 keV the maximal depth for such a line would be tcrsf ^ 0.23, 
assuming a line with a Gaussian optical depth. This is roughly a 
factor of two less than the depths of the lines seen for example 
in A 0535+2 69 dCaballero et ap |2007| or4U 1907+09 ( |Rivers| 
20 1 0) > . As 4U 1909+07 is clearly a young neutron star, 
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as proven through the strong pulsations and the early type stel- 
lar companion, it also possesses a strong magnetic field, which 
in theory must lead to the formation of CRSFs. As theoretical 
simulations show, however, it is possible that these features (es- 
pecially the fundamental line) can be filled up again by photon 
spawning and are therefore not detectable in the observed X-ray 
spectrum ( Schonher r et al.| |2007[ ). Additionally, the geometry of 
the magnetic field as well as the viewing angle onto the neutron 
star both can have strong effects regarding line shape and depth. 
Work is still in progress in understanding and quantizing these 
effects. The stringent upper limit on the CRSF depth is therefore 
very important to better understand the physics and magnetic 
fields of neutron stars. 
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